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ABSTRACT 

We update the limit from the 90 GHz PIQUE ground-based polarimeter on the magnitude of any 
polarized anisotropy of the cosmic microwave radiation. With a second year of data, we have now 
limited both Q and U on a ring of 1° radius. The window functions are broad: for E-mode polarization, 
the effective I is {is) ~ 191^132- We find that the E-mode signal can be no greater than 8.4 /xK (95% 
CL), assuming no B-mode polarization. Limits on a possible B-mode signal are also presented. 

cosmology: polarization — cosmology: 



Subject headings: cosmology: cosmic background radiation 
observations 



1. INTRODUCTION 

The current limits on cosmic microwave background 
(CMB) polarization restrict the amplitude of its 
fluctuations to less than 10 /iK at 95% CL. At 
large angular scales, Keating et al. (2001) limit the 
amplitude to 8 /iK. At subdegree angular scales, the 
constraint from Hedman et al. (2001, hereafter HOI) 
is 10 /iK, while at arcminute scales, Subrahmanyan et 
al. (2000) set a limit of 10 /tK. Present estimates of 
the peak polarized fluctuation amplitude are 6 /iK 
at an angular scale of 0?2 (£ ~ 950). These 
estimates are based on parameters gleaned from CMB 
temperature anisotropy measurements (e.g. Pryke et 
al. 2002; Jaffe et al. 2001; Wang, Tegmark, & Zaldarriaga 
2001). While CMB polarization has yet to be detected, 
its characterization will complement CMB temperature 
anisotropy data and impact our understanding of: 
gravitational waves from the inflationary epoch (e.g. 
Turner 1997, Caldwell, Kamionkowski, & Wadley 1999); 
peculiar velocities at the surface of last scattering 
(Zaldarriaga & Harari 1995); the nature of primordial 
perturbations (e.g. Spergel & Zaldarriaga 1997); 
primordial magnetic fields (Kosowsky & Loeb 1996); 
and cosmological parity violation (Lue, Wang, & 
Kamionkowski 1999). Here, we report improved limits 
derived from new data from the 2001 observing season of 
the Princeton IQU Experiment (PIQUE) at 90 GHz. We 
combine the new data with data from the first observing 
season, and also present a reanalysis of those earlier 
data. These data pass extensive checks for systematic 
contamination. Future publications will report results 
from a 40 GHz polarimeter also deployed during the 2001 
observing season, and give details of the instrument. 

2. INSTRUMENT, OBSERVATIONS AND CALIBRATION 

PIQUE has been described previously (HOI). The 
results reported here are from PIQUE's broadband 
90 GHz correlation polarimeter, which underilluminates 
a 1.2 m off-axis parabola (WoUack et al. 1997), resulting 
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in a beamsize of 0?235. The 84-100 GHz bandpass is 
divided into three subbands called SO, SI, and S2 (HOI). 
Observations are made of a ring of radius 1° around the 
NCP; the telescope site is Princeton, NJ. 

The polarimeter observed the sky from 2000 January 
19 to 2000 April 2 and from 2000 December 19 to 2001 
February 28. These two observing seasons yielded 810 hrs 
and 660 hrs of raw data, respectively. 

The scanning strategy is designed to permit null tests for 
checking sensitivity to systematics. During both observing 
seasons, the telescope alternated between two azimuth 
positions at fixed elevation. Data from the two positions 
are differenced to remove sensitivity to DC offsets. For the 
first season these azimuth positions were ±0?93, and the 
elevation was 41?0. The telescope therefore measured 
(as defined by the lAU) for two regions separated by six 
hours in right ascension (RA) on the ring of declination 
89°. For the second season the azimuth positions were 
±1?31, with elevation 40? 3, so the telescope measured +U 
for two regions separated by twelve hours in RA on the 
same ring. The azimuth chop period was 13 seconds until 
2001 January 23, at which point it was doubled. 

The polarimetry channels are calibrated to 10% using 
a nutating aluminum flat (HOI, Staggs et al. 2002). 
Constant elevation scans of Jupiter are used to determine 
pointing accuracy and map the beams. For the second 
observing season the measured beam FWHM are 0?235(7) 
in co-elevation and 0?233(7) in elevation, in agreement 
with measurements from the first season. The absolute 
pointing offsets in elevation are smaller than 0?03. 
However, early in the second observing season, the 
encoder suffered a misalignment so that the azimuth offset 
increased from -t-0?03 to -f0?08. (No concomitant change 
in beamshape was observed.) Note that we are able to 
neglect this in the analysis because the overlap between 
the ideal and misaligned beams is still 85%. In fact, 
simulations indicate the misalignment has less than a 2% 
effect on our derived limit. 

3. DATA REDUCTION 

The 1470 hours of data from the two observing seasons 
include: 383 hours of data taken while the telescope was 
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Fig. 1. Binned data are shown in tliermodynamic 
units for each observation season. Tlie three frequency 
channels have been coadded and the one sigma errors include 
small correlations among the channels. For clarity only 24 
approximately beam-sized bins are shown: the actual analysis 
uses 144 independent bins. Due to PIQUE's differencing 
strategies, a 1 Jy point source with 20% polarization located 
at a = 12'', 3 = 89° generates the point source responses 
shown in the two panels. No such point source is known, and 
Toffolatti et al. (1998) estimate there are fewer than 200 point 
sources brighter than 1 Jy in the entire sky at 90 GHz. 



slewing between the desired scan positions, 339 hours of 
data corrupted by known electromechanical failures, and 
59 hours of data in isolated fragments less than 12 hours in 
length. The remaining 685 hours of data contain 186 hours 
of data corrupted by meteorological phenomena (clouds), 
which are identified using the selection criteria described 
below. 

As discussed in HOI, the correlation polarimeter suffers 

a small sensitivity (< 0.5%) to total power because signals 
reflected from the input of the amplifier in one arm 
can couple into the other arm through the orthomode 
transducer. Therefore, the distributions of correlation 
coefficients 

< SiSj > 



< Sf >< S] > 

between pairs of polarimetry channels (i,j) display 
large positive tails due to periods of rapid atmospheric 
fluctuations. Data corrupted by clouds are removed 
by requiring the coefficients to be less than certain 
thresholds. Such selection criteria are determined based 
upon a data set designed to be insensitive to real 
astronomical polarized signals: the quadrature data (for 
which data from each scan position are split into halves 
and differenced). The Cjv are generated as averages over TV 
chops; varying N varies the time scale probed. For purely 
Gaussian noise, the shapes of the resulting distributions 
of Cjv for the whole data set depend only on N. In order 
to avoid using additional cutting measures (such as the 6- 
hr null test used in HOI), the CAr-selcction technique has 
been refined from that used for HOI. Here, two time scales 
are used rather than one. First we calculate the average 
correlation coefficients for segments of the time series 40- 
70 minutes long (specifically, N = 200). Segments with 



coefficients larger than a threshold of 0.20 arc removed. 
Next, coefficients are calculated for N = 50, for which the 
cut threshold is 0.32. These thresholds are at 2.5cr and 2a, 
and are chosen so that either cut alone removes ~ 20% of 
the data. The combined cuts remove 27% of the data. 
Null test results are not sensitive to the exact values for 
the thresholds. 

The 307(192) hours of data surviving the cuts detailed 
above from the first (second) observing season are parsed 
into 144 bins based on the Local Sidereal Time (LST) 
when the data were taken, following the same procedures 
outlined in HOI. These data are plotted in Figure 1. 
Offsets on the order of a few hundred /iK are removed from 
each polarimetry channel for each "deployment" (a period 
of > 12 hours bracketed by periods when the instrument 
was tarped). The results are not sensitive to the exact 
number of offsets removed. Table 1 presents the results 
of the null tests described in HOI, using the new selection 
criteria. The distribution of these null tests is consistent 
with noise, demonstrating that the data do not suffer from 
residual atmospheric contamination. 

4. DATA ANALYSIS 

The likelihood of a model given a data vector x 
is £ oc exp(— x'^C~-'x/2)/|C|'/^, where the covariance 
matrix C sums both theoretical correlations from the 
model (signal) and correlations specific to the experiment 
(noise): C = S -|- N. The analysis of the combined 
2000-2001 data uses x with 864 elements comprising 144 
spatial pixels measured in three different frequency bands 

Table 1 

Results of consistency tests. 



Year 


Test 


SO" 


SI" 


S2" 


2000 


Quad** 
Hl-H2'^ 


0.96 

0.18 


0.25 

0.14 


0.19 

0.33 


Pattern'* 
Si-Sj« 


0.46 
0.62 


0.26 
0.83 


0.14 
0.83 


2001 


Quad'' 
Hl-H2= 


0.38 
0.17 


0.61 
0.61 


0.28 
0.58 


Pattern'* 
Si-Sj« 


0.80 
0.47 


0.12 
0.66 


0.38 
0.04 



"Each numerical entry gives the probability of exceeding 
the for the given frequency channel. 

''The quadrature test uses data from each scan position 
(east and west) split into two halves and differenced to yield 
the quantity (Ei - E2)/2 + {Wi - W2)/2. 

'^Data from the second half of each season are subtracted 
from the first half. 

** Pattern nulls are generalizations of the 6- hour null test 
from HOI, and are data sets constructed from the various 
combinations of the data that should be zero given the 
differencing scan strategy. If dt is the measured signal 
at Local Sidereal Time t in hours, then for 2000 these 
combinations are dt — dt+e + dt+12 — c^t+is, while for 2001 
these combinations are dt + dt+i2- 



'^Data diff'crcnced between two channels, 
entries are SO-Sl, S0-S2 and S1-S2. 
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over the course of two observing seasons. The noise 
matrix N encodes the variances for each pixel for each 
channel and also accounts for interchannel correlations 
from both atmospheric fluctuations weakly coupled into 
the polarimctcr channels and correlated gain fluctuations 
in the cryogenic amplifiers. The interchannel correlation 
coeflicients are < 8% on average, and smaller for the two 
most sensitive channels. (S2 has just 20% of the total 
weight.) The noise matrix docs not include pixel-pixel 
correlations, since no such correlations arc observed in 
the time series data. An independent analysis using a 
288-element data vector, for which data from the three 
frequency channels are combined for each spatial pixel, 
with errors calculated to account for the interchannel 
correlations, yields consistent results. 

Since PIQUE measured Q in 2000 and U in 2001, the 
signal covariance matrix takes the form: 



^^({QQ) {QU)\ 

\{QU) {UU)J 



(2) 



where (QiQj) (or (UiUj)) represents the theoretical 
correlation between two spatial pixels from the 2000 (or 
2001) data set, and (QiUj) encodes correlations between 
pixels from different years. Note that the Q are sums of Q 
separated by six hours in RA and the U are differences of U 
separated by 12 hours in RA. The expression for {QiQj) is 
given in HOI. Following Zaldarriaga (1998), the expression 
for (QiUj) is given in terms of the E- and B-mode angular 
power spectra Cf and Cf by 



{QiUj) = E 



tm 



[Cf +Cf]W^2A<t>i3), (3) 



where (j)ij = cos ^{fii-nj) is the lag. The window function 
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Fig. 2. — Zero-lag window functions for El- and B-modes (solid 
and dashed lines) are shown in the middle panel for the 2000 

observation season and in the lower panel for the 2001 season. 
The top panel shows the limit on Te assuming Tb = for 
the combined W-band observations. For comparison we also 
plot (dashed line) E-mode predictions from the best-fit model 
in Pryke et al. (2002), as well as B-mode predictions (shown 
in the top panel multiplied by ten) assuming the same model 
with T/S = 1. 



Wi2^e has the form: 

Wi2,e = ^(Sl^^)'-Fi,£mf2,^m cos(m<^y ), (4) 

m 

where the -F{i.2}jm ^-re given in terms of associated 
Legendre polynomials evaluated at the ring radius = 1°. 
Here, the beam function is: 

^^em) -^siiH^ 2 4 J (m7r/144)2 ^ 

(5) 

where a = 0?10 for the PIQUE beams. Similarly, the 

expression for {UiUj) is: 

^^^[CfW2A^,j)+CfW,Aci>,j)], (6) 



{UiUj) = 



where the Wi^e and W2^e are the associated window 
functions given in HOI. In this case, the beam function 
is given by 



(S, 



UU\2 
hn ) 



/RQ;7^,2sin(m7r/2) 



(7) 



cos(m7r/4) 

This formalism allows for combination of the two data sets; 
however, separate analyses of the Q data and the U data 
are also presented. We plot the zero-lag window functions 
in Figure 2. 

The likelihood analysis proceeds by considering flat 
angular spectra, such that l{l + 1)0-^/211 = T^, where 
X = E,B. Since the amphtude of is predicted to 
be much smaller than that of Cf, we first find the limit 
Te on Te under the assumption Tb is identically zero"'^. 
This limit is compared to predictions in Figure 2. Next, 
joint upper limits {Te,Tb) are determined by finding the 
constant contour of £{Te, Tb) enclosing 95% of the volume 
of £. This is repeated for each year separately and for the 
combination of the two years. The 95% confidence level 
(CL) upper limits are shown in Table 2 and the normalized 
likelihoods are shown in Figure 3. The null data sets are 

1 The limit Te is found by integrating £(Tb,0); the result is RJ 30% 
higher if C{T^,0) is integrated. 

Table 2 

^}')% CoxFiDEKCE Limits 



Year 


Data* 


Te^ 

(mK) 


Te" 

(mK) 


Tb" 

(mK) 


2000 


CMB 
Quad 
(Hl-H2)/2 


12.7 
10.3 
14.9 


15.8 
13.9 
19.7 


14.7 
12.9 
17.4 


2001 


CMB 

Quad 
(111-112)/ 2 


10.4 

10.1 
11.2 


15.9 

16.2 
17.2 


17.8 

18.0 
liJ.l 


2000 

-1- 
2001 


CMB 
Quad 
{Hl-H2)/2 


8.4 
6.1 
8.5 


11.2 
8.3 
11.5 


11.5 
8.4 
11.6 



''The two null tests are described in Table 1. 
'•The Umit Tg is found assuming Tg = 0. 

'^The limits Te and Tb are determined simultaneously 
by finding the contour of constant likelihood enclosing 95% 
of the volume. 
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treated in an identical manner; results are tabulated in 
Table 2 and plotted in Figure 3. The limits in Table 2 
do not include calibration errors. Note that for PIQUE's 
broad window functions, the 5% beam errors only add 2% 
errors in quadrature with the 10% calibration errors. 

5. DISCUSSION 

The main result here is a new constraint on the amount 
of polarized anisotropy in the CMB at sub-degree angular 
scales. The result derives from combining data on Q 
from our first campaign with U from our second; in so 
doing, important information on the Q-U cross-correlation 
is included. We have summarized this result as a 95% 
CL limit of 8.4 fjK on E-modes. Given PIQUE's window 
functions and current theoretical predictions (Figure 3), 
wc might expect a signal of a few /iK. The likelihood 
for our CMB data (Figure 3, bottom panel) is consistent 
with this expectation. When wc fit to an offset lognormal 
distribution (Bond, Jaffc, & Knox 2000), wc find a central 
value of 2.2 /iK , a variance of (16 /iK )^ and a noise- 
related offset of 18 /iK^. 

The results presented here have been checked with 
two independent likelihood analyses and supported by 
extensive simulations. Our new selection criteria, as 
described above: are better able to deal with instrumental 
effects on different time scales; allow us to discard the 6-hr 
null test criterion we previously used; and work for both 
data sets together. However, applying the new criteria 
just to year 2000 data, we find a weaker 95% CL limit on 
polarized CMB anisotropy than in HOI: 12.6 /xK rather 
than 10.3 /xK. 

We have performed a variety of simiilations to address 
the probability of such a change. Recall that our result for 
the year 2000 data was essentially unchanged by relaxing 
the null test cut and allowing in 80 extra hours of data, for 
a total of 330 hours. For these simulations, we start with 
roughly a 330 hour data set which is pure noise, generated 
assuming the actual weights for each period of the data 
set. We then investigate how cutting the data can change 
the derived limit. From this we find that although a) the 
probability of the limit not worsening when 330 hours is 
reduced to 250 hours is only about 1%, b) a 10 /xK error 
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Fig. 3. — Normalized likclifioods versus fiat-bandpower Te 
under the assumption Tb = 0, as described in the text, for 
each year of data, and for the combined data. In each panel, 
likelihoods for the CMB data are plotted along with likelihoods 
for two of the null data sets described in Table I. 



from 330 hours of data is within one standard deviation of 
what is expected from pure noise (given our experimental 
weights); and c) the expected change in the limit from 
removing chunks of the data to get to 300 hours of data, 
(our final sample for year 2000), is about +1.1 /iK with 
a standard deviation of about 1.6 /xK. Thus our observed 
change is again within one standard deviation of what is 
expected. 

We thus conclude that fluctuations alone can account 
for the change in the limit derived from the year 2000 data 
under different selection criteria. This is so even though 
the former cut on the null test, which we have shown is 
not needed, was probably too restrictive. 

The expected signal, even including foregrounds, is 
still smaller than our new limit. Given the multipole 
range probed by PIQUE, this result provides the tightest 
constraint yet on the polariation spectrum predicted from 
primordial density fluctuations. 
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